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USE OF HEXAFLUOROETHANE TO SIMULATE THE INVISCID 
REAL-GAS EFFECTS ON BLUNT ENTRY VEHICLES 

By James L. Hunt, Robert A. Jones, and Kathryn A. Smith 
Langley Research Center 

SUMMARY 

Results are presented of an experimental and theoretical investigation of the use of 
equilibrium hexafluoroethane (C2F0) as an undissociated test gas in a wind tunnel to sim- 
ulate the inviscid aerodynamic characteristics of blunt bodies in high-speed flight where 
dissociation occurs. The results indicate that the use of C2F0 as a test gas in wind tun- 
nels is a practical and relatively simple and inexpensive method of obtaining test conditions 
with large normal- shock density ratios which is the necessary condition for blunt-body 
inviscid similitude. Equations for the thermodynamic and transport properties of C2F0 
are also included. 

INTRODUCTION 

A study was made in reference 1 of the simulation of reentry conditions which can be 
obtained in ground test facilities. This study concluded that the high Mach number, high- 
enthalpy characteristics of reentry conditions were not duplicated in operational ground 
facilities. These high-enthalpy conditions result in highly compressible real-gas flows. 

As a blunt body enters a planetary atmosphere at hypersonic speeds, the gas molecules 
that pass through the bow shock are excited to higher vibrational and chemical energy 
modes. This lowers the specific-heat ratio of the gas below the free- stream value, if it 
is assumed that equilibrium exists and provided that dissociation is not driven near com- 
pletion (both of these conditions aie usually met for blunt-body entry). A large amount of 
the energy that would have gone into increasing the static temperature behind the bow 
shock is thus used in exciting vibration energy levels or dissociating the gas molecules. 

As additional energy is absorbed by the gas molecules entering the shock layer, the con- 
servation laws and the thermophysics of the gas dictate certain changes in the forebody 
flow. The static temperature, speed of sound, and velocity in the real-gas shock layer are 
reduced. The static pressure is increased slightly. The density is increased consider- 
ably, and the shock-layer thickness is reduced in proportion to this increase. 




Use of a substitute gas in simulating the inviscid aerodynamics on blunt and slender 
bodies at moderate to large angles of attack in these highly compressible real-gas flows 
without simulating enthalpy was presented in reference 2 for blunt bodies and applied to 
slender bodies at high angles of attack in references 3 and 4. The contentions therein, 
references 2, 3, and 4, were as follows: (1) For hypersonic flow conditions, aerodynamic 
characteristics of blunt bodies are determined primarily by the inviscid flow field which, 
in turn, is most influenced by the shock density ratio (an indicator of the amount of real- 
gas effects present) in the blunt, forebody regions where subsonic flows persist and by the 
local isentropic exponent and local Mach number aft of the forebody in supersonic regions. 
(2) For ground tests at hypersonic Mach numbers where no dissociation occurs, the density 
ratio across a shock is determined by the specific -heat ratio of the gas; thus, ground-test 
data obtained at hypersonic Mach numbers at various values of specific-heat ratio (or 
shock density ratio) could be used to simulate the trends of equilibrium real-gas effects 
expected in flight. 

Due to dissociation of the gas in the shock layer, shock density ratios as high as 15 
and 20 are encountered during planetary entry; however, wind tunnels using air or nitrogen 
and operating at perfect gas conditions are limited to shock density ratios of less than 6. 

By using tetrafluoromethane |CF 4 ) as the test gas in a Mach 6 tunnel (ref. 2), density 
ratios as high as 12 have been obtained at tc f al temperatures less than 800 K with no dis- 
sociation. These high density ratios are possible because of the very low ratio of specific 
heats of the large CF 4 molecules; this CF 4 tunnel operates virtually as if it had a perfect 
gas flow with a specific-heat ratio of 1.12. 

The high cost of tetrafluoromethane (CF 4 , $8.80 per lb, July 1974) and the need to 
extend the density- ratio- simulation range above 12 at hypersonic Mach numbers has led 
to an exploratory experimental and theoretical investigation of the possibilities of using 
hexafluoroethane ^Fg, $3.00 per lb, July 1974) as an undissociated test gas in a wind tun- 
nel to simulate the aerodynamic characteristics of blunt bodies. The purpose of this 
investigation is, therefore, to determine whether or not the use of C 2 F 0 as a test gas in 
wind tunnels is a feasible method of obtaining test conditions with large normal-shock den- 
sity ratios (the necessary condition for blunt-body inviscid similitude). 

SYMBOLS 

A,B,C,D,E,F constants in saturation pressure equation (eq. (A15)) 

a,b,c,d,f constants in equation for specific heat at zero pressure (eq. (A3)) 

a a speed of sound 
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specific heat at constant pressure 


specific heat at constant volume 


ideal-gas specific heat at constant volume 


enthalpy 


thermal conductivity 


thermal conductivity at atmospheric pressure 


Mach number 


molecular weight 


free-stream Reynolds number (based on base diameter) 


pressure 


gas constant (eq. (Al)) 


radius 


entropy 


temperature 


velocity 


specific volume 


compressibility factor, pV/RT 


constants in equation of state (eq. (Al)) 


ratio of specific heats, Cp/Cy 


effective value of y 
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isentropic exponent (eq. (A12)) 
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A 

shock standoff distance 


e 

inverse density ratio, 


M 

viscosity 


M* 

viscosity at atmospheric pressure 
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viscosity parameter, 


P 

density, 1/V 


Pr 

reduced density, V c /V 


0 

proportionality constant (eq. (A17)) 


n 

constant in equation (A18) 


Subscripts: 



c 

critical point 


r 

reference conditions 


s 

saturation 

C 

0 

free-stream stagnation conditions 


1 

free- stream static conditions 
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static condition behind normal shock 
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stagnation condition behind normal shock 
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SIMULATION OF REAL-GAS EFFECTS ON BLUNT BODIES 

At supersonic and hypersonic speeds, the aerodynamic characteristics, including 
the lift, drag, moments, and stability, of blunt (high drag, low fineness ratio) configurations 
are determined almost exclusively by the forebody flow field. Several previous investi- 
gations (refs. 5 to 8) have shown, for example, that the drag of such configurations in a 
constant y flow (air) is insensitive to Mach number. Other investigations (refs. 2, 7, 8, 
and 9) have shown that for very high speeds such configurations have aerodynamic charac- 
teristics that, although essentially independent of Mach number, are strongly dependent on 
real- gas effects such as excitation of vibrational energy levels or dissociation. These 
real-gas effects have been shown to correlate as a function of density ratio across the 
strong bow shock ( P\/P 2) * or klunt configurations. 

The changes in aerodynamic characteristics due to real-gas effects associated with 
high-speed flight (characterized by large shock density ratios) are primarily the results 
of changes in surface pressures acting on the forebody. The surface pressures are 
affected by a change in shock density ratio (real- gas effects) in two ways: first, the level 
of pressure at the stagnation point is changed, and, second, distribution of surface pres- 
sure relative to stagnation- point pressure is changed (for details, see refs. 2 and 10). 

The investigation of reference 2 was restricted to bodies on which sonic velocities occur 
at the point of maximum body radius. This is not to say that the density ratio is not a 
simulation parameter on blunt bodies where the sonic line has moved forward of the point 
of maximum body radius. 

One way to simulate inviscid real- gas effects on blunt bodies is to test in a substi- 
tute gas flow which provides the correct value of density ratio without any dissociation. 

The shock density ratio in the substitute gas is a function of both Mj and y. For an 
ideal gas (ref. 11), the expression for the density ratio across a normal shock is 


p 2 (y + l)Mj 2 
P 1 (y - l)M t 2 + 2 


( 1 ) 


An effective value y e ff can be defined as the ideal or substitute gas value which gives 
the correct normal-shock density ratio when the free-stream Mach number in equation (1) 
is used. The distribution of density ratio with shock angle in flight is such that matching 
the normal-shock density ratio also insures matching of the oblique-shock density ratio 
for shock angles above 50°. At large hypersonic Mach numbers, the density ratio becomes 
a function of y only; however, as y decreases, the Mach number at which this condi- 
tion is reached increases. 
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CHARACTERISTICS OF HEXAFLUOROETHANE 

Hexafluoroethane (Freon- 116, C 2 Fg) is a nonlinear polyatomic molecule. Its 
17 vibrational energy modes (ref. 12) display anharmonic oscillator characteristics hav- 
ing 11 fundamental frequencies, 6 of which are twofold degenerate. The acute tempera- 
ture dependence of the vibrational energy of this molecule reveals itself in the specific 
heats at constant pressure and volume (c p and c v ) and their extremely low ratio y 
along with the isentropic expansion exponent (y ie j. 

Gaseous undissociated hexafluoroethane is colorless, odorless, and nontoxic. Toxic 
decomposition products may be formed upon thermal or electrical breakdown. However, 
C 2 F6 is among the most stable of organic compounds (refs. 13 and 14). Reaction with 
quartz does not occur below about 873 K, and homogeneous decomposition is not appre- 
ciable at 1115 K (ref. 13). Some of the physical properties of C2Fg (obtained from 


E. I. du Pont de Nemours & Co.) are: 

Molecular weight 138 

Boiling point at 1 atm (101.3 kN/m^), K 194.8 

Freezing point, K 172.4 

Critical temperature, K 292.7 

Critical pressure, atm 29.4 

Critical volume, cm^/mol 225 

Critical density, kg/m^ 611 

Density, liquid at 200 K, kg/m3 1587 

Density, saturated vapors at boiling point, kg/m3 9.01 

Specific heat, liquid at 200 K, kJ/kg-K 0.971 

Estimated specific heat, vapors at 1 atm and 298 K kJ/kg-K 0.761 

Estimated specific-heat ratio at 1 atm and 298 K 1.085 

Heat of vaporization at boiling point, kJ/kg 117.04 

Vibrational relaxation time at 298 K and 1 atm (ref. 15), sec 2.7 x 10"® 

Thermal conductivity at 298 K and 1 atm, W/m-K 0.017 

Prandtl number at 298 K and 1 atm 0.63 

Viscosity at 298 K and 1 atm, N-s/m^ 1.4 x 10‘® 


PREDICTED AERODYNAMIC SIMULATION IN C2F 6 

} 

In order to duplicate large hypersonic -flight density ratios in a wind tunnel without 
dissociation, the test gas must have a low ratio of specific heats and a relatively high ? 

degree of thermal stability (will not dissociate over the operating temperature range). 

Also, nonequilibrium and condensation in the test flow are highly undesirable; therefore. 
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the test gas should have a vibrational relaxation time that is low compared with the time 
required for the flow to span the shock standoff distance and a vapor pressure that is 
above the static pressure of the test flow. On the basis of the above criteria, the physical 
characteristics (see section entitled "Characteristics") indicate that hexafluoroethane 
may be a promising candidate for high-density-ratio flow simulation. 

In order to evaluate the suitability of hexafluoroethane as a test gas for high-density- 
ratio simulation, a computer code was developed (using the equation presented in the 
appendix with the shock-crossing equations of ref. 2) to calculate equilibrium conditions 
for isentropic expansion and flow across normal and oblique shocks in C2Fg. The struc- 
ture of this code, along with the shock crossing mechanics, is identical to that given for 
a CF4 code in reference 2. Some of the flow variables of interest which w ire computed 
from the computer code have been plotted as functions of pertinent parameters for the 
ranges of stagnation chamber pressure and temperature studied for stream Mach numbers 
from 3 to 7. A few of the variables of interest are given for CF4 at the design stagna- 
tion condition of a CF4 wind tunnel for comparative purposes. The stagnation enthalpy 
Hq of C2Fg, which is practically independent of stagnation pressure in the range from 
1034 to 1724 N/cm^, is given as a function of stagnation temperature Tq in figure 1. 

Plots of u x , P! u^Pq, Tj/Tq, Pl /p 0 , y v Y 2 /Y v P 2 / P l’ p l u l/ p l> P 2 U 2/ P 2’ 

Pg/pQ, and P3/P1 against Mj for C 2 F6 are presented in figures 2 to 12 for Tq = 478, 
588, 700, and 811 K. Specific -heat ratio, the unit Reynolds number, and the Reynolds 
number behind a normal shock are given for stagnation pressures of 1034, 1378, and 
1724 N/cm^. The other flow variables are independent (within the accuracy of the plots) 
of stagnation pressure in the range from 1034 to 1724 N/cm^. 
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Figure 4 shows the static temperature dropping below the critical at Mach 3.8 for 
Tq = 478 K, at Mach 5.2 for Tq = 588 K, at Mach 6.4 for Tq = 700 K, and above Mach 7 
for Tq = 811 K. However, the pressure required for condensation (vapor pressure) is 
much higher than the static pressure (fig. 5) over the Mach number range and stagnation 
conditions of this study; therefore, condensation does not occur. 

The similitude principle adhered to herein is to simulate the inviscid aerodynamics 
on hypersonic blunt bodies by matching the density ratio across the normal shock. For a 
given Mach number the normal- shock density ratio is shown to depend only on the y 
of the gas media in equation (1). Figure 6 shows that for a temperature of r ’l K the y 
C2F0 remains below 1.10 up to a Mach number of 7, whereas it exceeds 1.20 for CF4 at 
these conditions. A comparison of figures 6(a), 6(b), and 6(c) shows that the free-stream 
y becomes independent of pressure when the stream Mach number exceeds 3. Calcula- 
tions indicate that the stream isentropic exponent ^y ig j becomes essentially equal to the 
specific-heat ratio (y) once the stream Mach number exceeds 4. The shock specific-heat 
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\ ratio V 2 /y\ is given in figure 7 as 0.985 at a Mach number of 6 for a total temperature 

of 811 K. This corresponds (fig. 6) to a free-stream y of 1.069 and a y behind the 
shock of 1.052. Figure 8 shows that for a temperature of 811 K a shock density ratio of 
17,6 is obtained at a Mach number of 6 (y g £ f = 1.059). In CF^ a density ratio of 12 is 
obtained at these same test conditions (y e ff = 1.12). Scanning the stagnation condition of 
figure 8 for a Mach number range of 5 to 7 gives a density-ratio range of from 13 to 20. 

Inclusion of viscous regions in the aerodynamic similitude of blunt-body flows 
(incompressible flow) requires matching local Reynolds number. The unit Reynolds num- 
bers for the previously stated range of test conditions are given in figure 9. They are 
much lower than that for air at the same stagnation conditions. The Reynolds number 
range for a C2F0 Mach 6 tunnel based on a model length of 0.3 m would be 1.0 x 10^ 
to 4.0 x 10®. 

The Reynolds number behind a normal shock in C2F5 is given as a function of Mach 
number in figure 10. The Reynolds number behind the normal shock in C2F0 at Mj = 6 
is lower (5 x 10®/m versus 4 x lO^/m) than that in CF4 at Mj = 6 for the same stagna- 
tion conditions (Tq = 478 K; p Q = 1724 N/cm 2 ), but of the same order of a Mach 6 air wind 
tunnel (5 x 10Vm) at normal operating conditions. The large Reynolds number behind the 
normal shock in CF4 and C2F6 is due primarily to the large density behind the shock. In 
flight, where real-gas effects are prominent, large Reynolds numbers also occur behind 
shocks because of the large density. Noting that in any simulation it is the local conditions 
which should be simulated and not necessarily the free-stream ones, C2Fg and other high- 
density- ratio (low y) substitute gases may come nearer to simulating the local Reynolds 
number on a vehicle at flight conditions where real-gas effects are prominent than high 
Reynolds number air or helium tunnels (ideal gas-high y). This is especially true if the 
local Reynolds number and local Mach number simulation capability is examined behind 
20° to 50° oblique shocks when the no? mal- shock- density- ratio conditions are matched, 
since the low y gases also give higher Mach numbers behind oblique shocks than do the 
high y gases. 

The shock density ratio and Reynolds number simulation capability of CF4 and C2F0 
(assuming 10.15-cm-diameter model and a realistic pressure range) is shown in figure 13 
and that of several facilities at the Langley Research Center, which are thought to be typ- 
ical of most facility capabilities, is shown in figure 14. Also shown in both figures arc 
two entry trajectories typical of very blunt high-drag vehicles. Both these trajectories 
are based on the assumption of chemical equilibrium in the shock layer. The Earth entry 
trajectory is for an angle of 90°, an entry velocity of 5600 m/sec, a ballistic coefficient of 
146 kg/m 2 , and a diameter of 0.67 m. The Mars trajectory is for an assumed atmosphere 
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of 2000 N/m2 (model 2, ref. 16), an entry angle of 90°, an entry velocity of 8700 m/sec, 
a ballistic coefficient of 23.5 kg/m2 f and a diameter of 5 m. The "knee" of the trajec- 
tories corresponds roughly to conditions of poak heating and peak dynamic pressure. The 
Reynolds number behind the shock would be a better simulation parameter than that of the 
free stream to use in figures 13 and 14; this would show both CF4 and C2F5 in an even 
more favorable light concerning the simulation of the two trajectories. Since the aero- 
dynamics of blunt bodies depends essentially only on *he shock density r ,’e aerody- 
namics of a blunt Mars entry vehicle could be simulated over almost entire *&rs entry 
trajectory (fig. 13) in C2Fg. An examination of figures 13 and 14 shows that wit*. C^Fg 
the ability to simulate shock density ratios which occur with high-speed atmospheric entry 
is increased from 14 to 20. 

EXPERIMENTAL TEST AND RESULTS 
Apparatus 

Exploratory experimental tests were conducted in the Langley pilot CF4 Mach 6 
tunnel (ref. 2) to evaluate the suitability of C£Fg as a test gas. The tunnel is an inter- 
mittent type that is supplied by two large bag-type accumulators. These accumulator s 
maintain a constant tunnel stagnation pressure for times as long as 3 min. The gas 
passes from the accumulators through a lead-bath heat exchanger (stainless steel tubes 
immersed in a bath of molten lead) into the stagnation chamber and through the tunnel 
nozzle. The tunnel's contoured axisymmetric nozzle was designed for a CF4 flow at 
Mach 6 for stagnation conditions of 811 K and 1724 N/cm2, The nozzle is 67 cm long with 
a 16.4-cm-diameter test section. The expansion (area) ratio of this nozzle is 2569. The 
tunnel exhausted into a 1200- m3 vacuum sphere which could be pumped to pressures as 
low as 25 N/m2. No attempts were made to reclaim the C2Fg, although it should be pos- 
sible to do so easily since it is a heavy gas with a relatively high critical temperature 
(292.7 K) and a relatively low critical pressure (29.4 atm). To keep the gas above the 
critical temperature during pumping from the storage bottles to the accumulators (to avoid 
consensation at the high pressures), the storage bottles ^nd transport lines were wrapped 
with strip heaters and heated to 316 K. This facility description and operational proce- 
dure verifies the practicality of using C2Fg as a test gas in a hypersonic wind tunnel. 

Experimental Tests 

Tunnel survey .- A pitot rake with six orifices spanning 11.3 cm of the 16.4-cm- 
diameter test section was used to survey the C2F0 test stream for three stagnation condi- 
tions. These conditons were: 
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Survey test 

P 0 , N/cm 2 

T 0 , K 

1 

948 

520 

2 

1265 1 

581 

3 

1565 

603 


The results of the pitot surveys are shown in figure 15. The uniformity of the Mach num- 
ber distributions are acceptable over the 10.3-cm span, and the level Mi * 5.4 is consist- 
ent with that predicted from the real-gas ^Fg^ computer code for the area ratio of the 
nozzle. The tunnel -wall static pressure in the ,.est section was measured during each of 
these tests. The Mach number corresponding to the static pressure is also shown in 
figure 15 and is reasonably consistent with that corresponding to the total pressure 
measurements. 

Hemisphere tests .- Tests were conducted on a 1- in. -diameter hemisphere with a 
pressure orifice at its stagnation point at essentially the same three test conditions as 
that of the pitot surveys. Shadowgraphs of the bow shock (fig. 16) along with the pressure 
at the stagnation point of the model were obtained. Normal- shock density ratios in the 
test section were determined from the shock standoff distance measured on the hemisphere 
(fig. 16) by using the following relation (ref. 2): 

p 2 _ 0.76 

P 1 " A/r 

The normal shock density ratio obtained from the shock standoff distance along with the 
Mach number (M ca j c j calculated from the stagnation point pressures and the corresponding 
calculated density ratios (P 2 /Pi) ca j c * or eac h °* these tests are: 


Survey test 

p 0 , N/cm 2 

T 0 , K 


M C alc 

Ocalc 

4 

903 

505 

14 

5.40 

14.3 

5 

1250 

580 

14 

| 

5.40 

14.3 

6 

1500 

i 

589 

15 

5.35 

15.1 


The Mach numbers are consistent with those obtained with the survey rake and the calcu- 
lated density ratios (fig. 15). Tne density ratios obtained fron he shock standoff distance 
(eq. (2)) are in reasonable agreement with those obtained from isentropic expansion 
normal- shock calculations (fig. 8). 

The influence of the density ratio on the shock shape of a hemisphere is illustrated 
in figure 17 where sketches of shock shapes in air, CF 4 , and C 2 Fg are given at density 
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| ratios of 5.6, 12, and 16, respectively. The shock v/ave for the higher density ratio 

I (^2^6) is much closer to the body. Furthermore, since it lies closer to the body, the 

| expansion waves from the body and the expansion fan from the corner bends it more 

l sharply and thus the shock angle reaches its limiting value more rapidly than for the lower 

f density ratio shocks. This comparison of the shock shapes indicates that even though the 

limiting shock angle is 1.2° for air at Mach 8, 9.4 for CF^ at Mach 6.1, and 10.7° for 
i ^2^6 at Mach 5.4, a downstream distance of many diameters would be required before the 

\ shocks cross. 

CONCLUDING REMARKS 

An experimental and theoretical investigation of the use of hexaflaoroethane 
, (Freon- 116, C 2 Fg) as a test gas in a wind tunnel to simulate the inviscid aerodynamic 

characteristics of blunt bodies in high-speed flight where dissociation takes place in the 
shock layer indicates the following: 

1. The use of C 2 Fg as a test gas in hypersonic wind tunnels is shown to be a practi- 
cal and relatively inexpensive method of obtaining low specific -heat- ratio (y) flows with 
large normal- shock density ratios without dissociation. 

2. Calculations indicate that a shock-density-ratio range from 13 to 20 could be 
achieved with C 2 Fg for free-stream Mach numbers from 5 to 7 with stagnation tempera- 
tures below 811 K. 

3. Tests made in a small pilot C 2 Fg tunnel confirmed real- gas calculations for test 
conditions in C 2 F^ Values of normal-shock density ratio as high as 15 were obtained in 
the pilot C 2 Fg tunnel at a stagnation pressure of 1500 N/cm^, a stagnation temperature of 
589 K, and a free-stream Mach number of 5.35. 

4. Since the aerodynamics ol blunt bodies depends essentially only on the shock den- 
sity ratio, the aerodynamics of a blunt Mars entry vehicle could be simulated over almost 
the entire entry trajectory in C 2 Fg. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 5, 1974. 
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APPENDIX 

THERMODYNAMIC AND TRANSPORT EQUATIONS OF 
GASEOUS HEXAFLUOROETHANE 

The thermodynamic property equations of hexafluoroethane ^Fg) are derived from 
the equation of state and the expression f ; the specific heat at zero pressure. These 
properties are essential in predicting isentropic flows, flows across normal and oblique 
shocks, and the inviscid flow fields over vehicles. The transport properties are essential 
when viscous phenomena such as boundary layers on a nozzle or vehicular walls are being 
considered. 

The equation of state of hexafluoroethane was obtained from E. I. du Pont de Nemours 
& Co. and is 


RT QtT -1 / 2 
V - 0 ’ V(V + fi) 


(Al) 


where p is given in N/m^, V in m3 /kg, and T in K and where 


R = 6.024225 x 10 JAg-K 

0 = 5.133877 x 10 -4 m 3 /kg 

a = 7.657512 x 10" 3 J-m3-K^/^kg2 


The specific heat at constant volume is 


i 
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) 


c v 




dV 


(A2) 


and Cy is the specific heat at zero pressure. A third-degree expression in temperature 
for Cy (obtained from E. I. du Pont de Nemours & Co.) is: 


Cy = a + bT + cT 3 + dT 3 



(A3) 
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APPENDIX - Continued 


where ci, is in J/kg-K and where 


a = 9.952795 x 10 JAg-K 


b = 2.608545 JAg-K 2 
c = -1.926454 X 10" 3 JAg-K 3 
d = 1.857347 X 10“ 7 JAg-K 4 
f = -3.473992 x 10 5 J-KAg 

Evaluating the integral term of equation (A2) with the equation of state, equation (Al), gives 


Cy — Cy + 


+ 30X^2 Jy_±I 

4(3 \ V . 


The entropy S expressed as a function of T and V is 


r v /v 

V i 9T, 


^ d T + \ 

J T r T J V 


where the subscript r represents reference conditions taken herein as 


S r = 443.504 JAg-K 
T r = 172.222 K 
p r = 2.468323 X 10 4 N/m 2 
V r = 0.41489 mVkg 

From the integration of equation (A5) the expression for entropy becomes 


S - S r = a In 


(-L) + b (T - T r ) ♦ |(x2 - T t 2) ♦ f(T3 . „*) - - 4j) 

J-LlA) - ± T r -3/2 + azlll J_5L_\ 

\v r - ?) 2/3 r \Vr +P) 2(1 \y + (i) 
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APPENDIX - Continued 




The enthalpy H expressed as a function of T and V is 



and becomes in integrated form 


H - H r . a(T - T r ) + S(t 2 - T r 2) + |(t 3 - T r 3 )t | (t 4 - T r 4 ) - t(l - Xj 

i( T ' Tr ) - v ^^” 1/2 ‘ V I/2 ) 

I— —I 

\v v r + (j; 


- |f (T-l/2 - T r -l/2) t 


RVj 

V r - /3' 


- 0RT 


y -0 v _ - /3 


i ^ - aT - 1/2 ^ 


3oT~ 1//2 
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(A8) 


where H r = 68152.7 J/kg. 

Other variables of interest are: 
Specific heat at constant pressure 


Cp = c v - 



Specific-heat ratio 



(A9) 


(A10) 
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APPENDIX - Continued 


Speed of sound 


- li e\ 1/2 - L v 2^ ] 1/2 = hh. l 2 . v 2 ^ V 

'Ws L \ 9V /Tj \_ C v (vvyj V 9 V /t[ 


(All) 


Isentropic exponent 


ie yd In pj g P [d P j s 


I T V 


(A12) 


where the partial derivatives in equations (A9) to (A12) are 


/ap\ _ R ( qiT~ 3//2 

\BT/ V V-/3 2V(V + /3) 


£P\ _ . RT + gT~ 1//2 (2V + /3) 

9V /T (V - /3) 2 V 2 (V + /3) 2 


(A1J») 


(A14) 


These equations are limited to temperatures above the saturated vapor line. The satura- 
tion pressure in N/m 2 as a function of temperature in kelvins (obtained from E. I. du Pont 
de Nemours & Co.) is 


log p s = A + | + C log T + DT + - 12 [i og (F - T) + 0.25527] + 0.52553 (A15) 


where 


A = 0.3762282521 X 10 2 
B = -1.031833 x 10 3 


C = -0.1297816132 x 10 2 
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APPENDIX - Continued 


D = 0.01321347 
E = -0.537560386 
F = 0.297317 x 10 3 


This equation is only applicable below the critical temperature (t c = 292.7 k). Also, 
equations (Al) to (A14) are limited to temperatures below those at which dissociation 
occurs (T < 1115 K). 


Transport Properties of Gaseous C2F0 

Viscosity .- An expression for the vapor viscosity of C2Fg at 1 atm of pressure was 
obtained from E. I. du Pont de Nemours & Co; the expression is 


M* = (l2.778\/f - 78.8) x 10“ 7 


(A16) 


where p* is in N-s/m 3 and T is in K. This equation has experimental data to support 
its validity over a temperature interval from 230 K to 420 K. The viscosity of C2Fg at a 
pressure of 1 atm determined from equation (A16) is plotted over a temperature interval of 
150 K to 800 K in figure 18. Even though the curve fit was only made for a temperature 
interval of 230 K to 420 K, equation (A16) behaves well over the extended range. To obtain 
the viscosity of C2F0 at pressures other than 1 atm, the residual viscosity is expressed 
as a function of the density, molecular weight, and critical constants of the substance 
(ref. 17) through a dimensional analysis treatment (ref. 15) resulting in the following 
expression: 


(p - ju*)| = 


r 1 / 6 



j 

R 


(A17) 


Experimental high-pressure viscosities for the gaseous and liquid phases available in the 
literature for 14 polar substances including five Freons ^CClgF, CCI2F2, CHC^F, CHCIF2, 
and C2CI3F2) were used (ref. 18) in conjunction with pressure-volume-temperature data 
to establish the constant a and the exponents i and j of equation (A17). Plots of 
(/i - p*)i j against (ref. 13) were found to be essentially the same for all 14 sub- 
stances; therefore, it was concluded that the exponent i of Z c in equation (A17) was 

1 zero. 
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APPENDIX - Concluded 

The analytical representation of the residual viscosity correlation for the 14 polar 
substances (ref. 18) which should apply to C2F0 as well are: 

<M - M*)$ = 16.56 x 10-8^1.111 fa < 0.10) (A18a) 

(M - j u*)£ = 0.607 x 10-8(9.045p R + 0.63) 1 - 739 (0.10 < p R £ 0.90) (A18b) 

-log[(p - #i*)*] = 100-6439-0.1005^-0 + 3 (o.9 < Pr < 2 .6) (A18c) 

where 0 = 0 for 0.9 < p„ < 2.2, 0 = 4.75 x 10~ 4 (p„ 3 - 10.65) 2 for 2.2 < p„ < 2.6, 

t1/6 , 9 

and £ = — = 0.023735 for CoFc and where p is in N-s/W. 

M‘/2 Pc 2/3 

Thermal conductivity .- An estimate for the thermal conductivity of C 2 Fg gas at a 
pressure of 1 atm was obtained from E. I. du Pont de Nemours & Co. as: 


k* = (9.591 x 1(T 5 )T - 0.01163 


(A19) 


where k* isinW/m-K. The estimate was made over a temperature range from 230 K 
to 420 K. The thermal conductivity at 1 atm of pressure is plotted over a temperature 
range of 100 K to 800 K in figure 19. 

The thermal conductivity relationship is extended to pressure other than 1 atm in 
the same manner as was the viscosity. Residual thermal conductivities k - k* for CF4 
were obtained from experimental data and correlated with the corresponding reduced den- 
sities (ref. 19) to produce the following expression: 

k - k* = 7.18 x 10-3|exp(p R ) - lj (A20) 

where k - k* is in W/m-K. Equation (A20) represents the thermal conductivity of CF4 
in the dense gaseous region and should also be fairly representative of the thermal con- 
ductivity of CgFg. 
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Figure 5.- Variation of ratio of free-stream static pressure to stagnation 
pressure with Mach number for C2Fg. (Independent of stagnation pres- 
sure from 1034 to 1724 N/cm2.) 
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Figure 14.- Density-ratio simulation capability of several facilities. 
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Figure 15.- Mach number distribution in pilot C 2 F 6 tunnel (axisymmetric nozzle). 












